Cooling a micro-mechanical resonator by quantum back-action from a noisy qubit 
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We study the role of qubit dephasing in cooling a mechanical resonator by quantum back-action. 
With a superconducting flux qubit as a specific example, we show that ground-state cooling of a 
mechanical resonator can only be realized if the qubit dephasing rate is sufficiently low. 
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A micro (nano) mechanical resonator (MR) can be 
cooled down through quantum back-action of coupled 
auxiliary mesoscopic systems. Environmental fluctua- 
tions induce both relaxation and dephasing processes to 
the auxiliary systems. Previous studies [E 0, d, 0, d, @, 0j 
[ill ] show that relaxation plays an essential role 
to dissipate the MR energy to the environment. In this 
paper, we address the function of environment-induced 
dephasing in the back-action cooling. 

Preparing quantum systems in their ground states is 
one way to initialize them for the implementation of 
many quantum protocols. If the MR is to be cooled 
to milli-Kclvin temperatures, novel cooling techniques 
other than dilution refrigeration are required. Experi- 
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mentally, cooling by an optical cavity (e.g. 
and references therein) and a superconducting single- 
electron transistor (SSET) has been demonstrated [15j |. 
but both techniques are still far from reaching the quan- 
tum ground state of the MR. Theoretical proposals have 
predicted the possibility of ground-state cooling of the 
MR 0, @, S, H i, H. The basic idea of these proposals 
resembles laser cooling of trapped ions. The dissipative 
auxiliary system (such as the internal levels of an ion, the 
optical cavity, the SSET, the quantum dot, etc.) acts as 
a structured bath of the system to be cooled (in our case, 
the MR). The relaxation and excitation rate of the aux- 
iliary system obeys the detailed balance relation deter- 
mined by the external bath temperature. Introducing a 
red-detuned drive modifies the dissipative nature of the 
auxiliary system: absorption processes in the auxiliary 
system and associated MR phonon emission processes 
dominate over the inverse processes, so that detailed bal- 
ance is broken. These phonon emission processes extract 
energy from the MR and dissipate it to the external bath 
through the subsequent spontaneous emission of the aux- 
iliary system energy quanta. Therefore, the limit of the 
cooling procedure is determined primarily by: (1) the en- 
vironmental noise of the auxiliary system and (2) the way 
the dissipative nature of the auxiliary system is modified 
by the drive. 



leads to both relaxation processes (i.e., energy-lowering 
or energy-rising processes) and dephasing processes. The 
cooling efficiency of trapped-ion or optomechanical cool- 
ing is only influenced by the relaxation process. In the 
case of cooling by driving a solid-state qubit @, 0, 13> EH] , 
pure dephasing can usually be neglected by biasing the 
qubit at (or close to) the degeneracy point where the pure 
dephasing rate is negligibly small. However, the first- 
order photon excitation by absorbing the energy from 
the low-frequency phonon and the linear drive vanishes 
at this point, only a small second-order photon interac- 
tion remains [tJ. Thus, the cooling becomes less efficient. 
This motivates us to bias the qubit away from the degen- 
eracy point and take the qubit dephasing into account. 
This consideration is especially important for solid-state 
qubits because the dephasing rate is much larger than the 
relaxation rate away from the degeneracy point and it is 
more difficult to suppress dephasing since it is associated 
with low-frequency noise. 

In this paper, we use the master-equation approach to 
investigate the ground-state cooling of the MR by a flux 
qubit. The influence of both qubit dephasing and relax- 
ation is studied. Assuming a 1// noise spectrum for the 
flux qubit, we show that ground-state cooling of the me- 
chanical resonator is possible under current experimental 
conditions. 

The system we consider is a micro-mechanical beam 
interacting with a superconducting flux qubit 17, [Til ], 
see Fig. [TJ In the x-y plane, a doubly-clamped micro- 
mechanical beam with an effective length Lq is incor- 
porated in a superconducting loop with three small- 
capacitance Josephson junctions. This mechanical beam 
can be created from an MBE-grown heterostructure 
coated with superconducting material or a self- 

supporting metallic air bridge. The fundamental vibra- 
tion mode of the beam can be well approximated as a 
harmonic resonator with oscillation frequency Ub- With 
a proper bias magnetic flux, two classically stable states 
of the 3-Josephson-junction loop carry persistent cur- 
rents in opposite directions. There is a finite tunneling 



The environmental noise acting on the auxiliary system rate A between the two classical persistent-current states 
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FIG. 1: (Color online) Schematic diagram of our setup. A 
doubly-clamped mechanical beam is incorporated in a super- 
conducting flux qubit (Josephson junctions are indicated by 
crosses). The initial bias in the loop is controlled by a mag- 
netic flux $ z o in the z-direction. A coupling magnetic field 
Bo is applied in the ^-direction, leading to a coupling of the 
motion of the beam and the qubit because the supercurrent 
leads to a Lorentz force on the beam. A microwave (MW) 
line introduces a microwave bias to the qubit loop. 



(throughout this article, we let ft = 1). This two-level 
subspace is far below the other energy levels and forms a 
flux qubit flc" " 



2l| . The qubit ground state \g) and excited 



state |e) are coherent superpositions of two persistent 
current states denoted by |0) (clockwise current state) 
and |1) (counter-clockwise current state). The energy 
spacing between the two eigenstates is ft — \J A 2 + e\ 
where £o = 2I p (<i> cx t — $0/2) is the energy spacing of 
the two classical current states (|0) and |1)), with <& cxt 
the external magnetic flux through the loop, I p the max- 
imum persistent current in the loop, and $0 = h/(2e) 
the flux quantum. A microwave line is placed close to 
the circuit and can be used as a microwave drive with 
frequency u>d acting on the flux qubit. The qubit Hamil- 
tonian is written as 

H q = (A/2)a x + {e /2)a z + Aa z cos(w d i) (1) 

where A characterizes the amplitude of the microwave 
drive and a x = \1) (0| + |0) (1|, a z = |0) (0| - |1) (1|. We 
assume that the drive is near resonant with the qubit 
\Sui\ / (ti>d +0) <C 1 (here, 5ui = ft — is the detuning 
between the qubit free energy and the drive), and that the 
drive amplitude satisfies A/ui^ <C 1. In the presence of a 
magnetic field B along the y-direction, the supercurrent 
generates a Lorentz force on the MR along the z- 
direction. This force couples the flux qubit with the mo- 
tion of the MR. The coupling Hamiltonian is g(a + a^)a x 
with g = BqI p LqSo- Here, <5o = \Jl/(2muJh) is the har- 
monic oscillator length (mean-square zero-point displace- 
ment). As shown previously using a semi-classical treat- 
ment fl(i ] , this configuration can serve as an "on-chip re- 
frigerator" for the MR under a proper drive: the Lorentz 
force produced by the flux qubit induces a passive back- 
action on the MR, damping its thermal motion. 

At low frequencies, the decoherence of the supercon- 
ducting flux qubit is dominated by 1// noise. This noise 



induces qubit relaxation and dephasing |22j, |23|, |24[ . In 



(excitation) rates 



rl -* satisfy the detailed balance 



relation /T^ — exp(f2 / k^Tp) , where T is the tem- 
perature of the external bath. Since £1 3> k^T in our 
case, we neglect in this discussion. The qubit pure 
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rate is almost proportional to the qubit 
3| in the vicinity of the degeneracy point so 



r(°) 



ae . 



(2) 



In the presence of a near-resonant microwave drive, in the 
rotating frame of frequency LOd, there are new eigenstates 
which are superpositions of the initial eigenstates |e) and 
I g). The new eigenstates are the qubit-microwave dressed 
states 25]. The dephasing of the undressed states also 
contributes to the relaxation and excitation of the dressed 
states. The relaxation (r^), the excitation (T^) and the 
dephasing rate (T v ) of the dressed states can be generally 
written as 

r i(T,v) = r i(T,v)R + r f.(T,v)D! ( 3 ) 
where r...R denotes the contribution from the relaxation 
process and r d from the dephasing process without 
drive. The derivation of the Lindblad-form master equa- 
tion of the driven qubit in the rotating frame leads to the 
following explicit expressions for the relaxation (excita- 
tion) rates: 



rf } (1 + cos <j)Y coth (J3 (top + o; d ) /2) ± I 
4 coth (/3fi/2) + 1 

cos cj)f coth (/3 (lo - Lu d ) 12) =p 1 
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4 coth (0CI/2) + 1 
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coth (/3Q/2) + 1 



= r<°>. 



<pT> = J-^'COS 



(0) 
(7) 

Equations (0])-([7]) are obtained by assuming a 1// noise 
spectrum [261 ] . Here, f3 = l/(k B T ), sin</> = A'/u> , 
coscf> = 5uj/ujq, ujq = V 8uj 2 + A' 2 , and A' — — Asin8 
with sin 9 = A/fl and cos 8 — eo/ft. Below the cutoff 
frequency uj c of the 1/ f spectrum, we assume the noise 
spectrum to be white noise. 

In the limit of weak qubit-resonator coupling, the Born 
approximation can be applied and one can trace out the 
qubit degree of freedom in the composite system master 
equation. This yields the following equation of motion 
for the phonon number 



-n(7m + 7q) + hmN + 7q Ay . 



(8) 



Its solution reads 



n (t) = nae-^J* + lmN + lqNq ( 1 - e^^^ 

7m + 7q V 



the absence of the microwave drive, the qubit relaxation 



(9) 
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and the stable state solution at t ^ l/(7m + 7q) is 



7/ 



7m + 7q 



(10) 



Here, 7 m is the damping rate and N — [exp (/3w b ) — 
the mean phonon number determined by the thermal 
bath, whereas 7 q = g[S zz (u>b) — S zz {— u b )] is the damp- 
ing rate and A q = S zz (-w b ) / (S zz (w b ) - S zz (-w b )) 
the mean phonon number determined by the qubit bath. 
Thus, the MR is effectively in contact with two baths: 
One is the real (external) thermal bath, and the other 
one is the structured bath formed by the dissipative qubit 
under drive (for simplicity, we call this the "qubit bath" 
in the following). Since the MR is coupled via the a z - 
component of the qubit, the MR decoherence related to 
the qubit bath is determined by the pair correlation func- 
tion S zz (lu) = J die Mt (a z (t) a z (0)), where the brackets 
denote the average over the steady state of the qubit 
that satisfies d{a XtV:Z ) / dt = 0. An explicit evaluation of 
S zz (oj) leads to 



S zz (w) 



cos 2 9 sin 2 i 



cos 2 9 sin 2 . 
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«+/ 2(uj + lu q ) +k, 
2k 



(11) 



Here, k = T^+Ti+2T V , k 4 



+ / 4 (uj — lu ) + K", 

r T +r b andd_ = r T -r x . 

Note that the condition to obtain this reduced master 
equation is the weak coupling between the MR and the 
two baths, i.e. j m + 7 q <C rj° ,u/ b . 

In the absence of the drive, 7 q = 0, and Eq. (JXOJ) shows 
that the final occupancy number of the MR is unchanged 
(n = N) . This is a natural consequence of the second law 
of thermodynamics. A slightly detuned drive breaks the 
thermal equilibrium and changes the final phonon num- 
ber of the MR. Figure [21 shows the dependence of the 
cooling efficiency n = N/n — 1 on the detuning Soj. It il- 
lustrates that cooling (77 > 0) is induced by a red-detuned 
(Su> > 0) drive: the red-detuned drive accelerates the 
photon emission process so that the MR energy is dis- 
sipated into the environment in the subsequent photon 
emission. In this way the mean phonon number of the 
MR is decreased. 

As an example, suppose we perform this cooling pro- 
cedure on a MR of length Lo = 5 um with oscillation fre- 
quency u>b — 10 MHz, spring constant k = 0.01 N/m, and 
quality factor Qm = 10 4 . The MR is coupled with a flux 
qubit loop by an in-plane magnetic field Bq — 6 mT. The 
qubit energy splitting is A = 5 GHz at the degeneracy 
point and I p = 400 nA. The coupling strength between 
the qubit and the MR is about 10.2 MHz. The initial 
environmental temperature is assumed to be Tq = 20 
mK. The qubit relaxation rate sw 4 MHz and the 

excitation rate 1^ ~ are almost independent of the 
qubit bias near the degeneracy point while the pure de- 
phasing rate depends linearly on the qubit bias in this 
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FIG. 2: (Color online) Inset: dependence of the final phonon 
number n on the scaled detuning 5u /u)b and the scaled drive 
amplitude A/uJb (with eo = 800 MHz). The plot range is 
limited from n = to n — 1 (i.e. the ground-state cooling 
regime). The region with n > 1 is shown in blue. Main panel: 
cooling power r\ = N/n—1 as a function of the scaled detuning 
Sio/ub between the microwave drive and qubit energy splitting 
along the dashed line in the inset. The other parameters are 
the same as those used to estimate the cooling limit in the 
text. 



regime so that » 0.008e [22], [23j. The cutoff fre- 
quency uj c is assumed to be 1 Hz. If the qubit is biased 
at £q = 1 GHz and driven by a microwave drive with 
A = 7.4 MHz, ui d = 5.056 GHz, the effective damping 
rate and the mean photon number of the qubit bath are 
7 q = 0.4 MHz, N q — 0.06 respectively. The steady state 
is reached after about 2.5 /is. The phonon number of the 
MR in the steady state is n = 0.12. Hence, ground-state 
cooling can be realized with this setup at the l//-noise 
level characterized by the parameters T i -* , rl°' ) , r^ -* , and 



T 



uj c given above which correspond to realistic values. The 
cooling power can be further improved by increasing the 
coupling magnetic field Bq. 

In order to achieve ground-state cooling, the power and 
detuning of the microwave drive should be optimized. As 
shown in the inset of Fig. [21 the final phonon number n 
(or the cooling power 77 = N/n — 1) does not depend 
monotonically on the drive magnitude. This is differ- 
ent from optomechanical cooling where a stronger drive 
leads to a stronger cooling effect. Qubit-assisted cooling 
requires a certain resonant condition and the qubit eigen- 
frequency in the rotating frame is modified by the drive. 
Previous experiments on cooling a MR by quantum back- 
action [15| show that the cooling limit is determined by 



the quantum fluctuations of the auxiliary system. In the 
following, we analyze how the qubit dephasing and relax- 
ation influences the cooling process in a different way. 

Since the decoherence of the qubit is sensitive to the 
qubit bias, the dependence of the final phonon number on 
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FIG. 3: (Color online) Final phonon number of the MR ver- 
sus scaled qubit bias. The three curves correspond to qubits 
with the same relaxation rate (r^ ~ 4 MHz) but different 
dephasing rates (from bottom to top a — (red), a — 0.008 
(blue), and a = 0.08 (black)). 

the qubit bias shown in Fig.[3]reveals the relationship be- 
tween cooling efficiency and qubit decoherence. Without 
consideration of dephasing (a = 0) , the cooling efficiency 
increases monotonically with the qubit bias. Taking a 
finite dephasing into account, there is an optimal point 
for the cooling efficiency as we increase the qubit bias. 
The cooling efficiency is decreased as the dephasing rate 
is increased. 

The behavior exhibited in Fig. [3] can be understood 
from Eq. (jlip . Increasing dephasing broadens the 
Lorentzians and decreases the difference between r-j- and 
The correlation spectrum S zz (cu) hence becomes 
more symmetric and the cooling effect is decreased. 
Therefore, N q increases with the increase of pure dephas- 
ing. Biasing the qubit to the degeneracy point helps to 
decrease qubit dephasing and hence decrease N q . How- 
ever, at the degeneracy point, 7 q = 0, photon excitation 
processes which absorb energy from the low-frequency 
resonator and the linear drive are not possible any more 
and the cooling cycle is stopped. Hence, there is an opti- 
mal qubit bias point that leads to a lowest final phonon 
number, see Fig. [3] 

The dressed-state relaxation and excitation rates have 
contributions from the qubit relaxation as well as pure 
dephasing. Physically, this is because the dressed states 
of the qubit in the presence of the microwave driving field 
are superpositions of the qubit eigenstates in the absence 
of the drive. The qubit pure dephasing contributes sig- 
nificantly to the decoherence process of the dressed qubit 
and hence modifies the cooling limit. The pure dephasing 
was not included in the previous studies of the trapped- 
ion cooling as well as the cooling of a MR by a quantum 
dot, optical cavity, or SSET. In these cases, the drive 
modifies the relaxation rate of the auxiliary system but 
leaves the order of the eigenlevels unchanged. The final 



relaxation process is not influenced by the pure dephas- 
ing. In the cooling schemes that use a driven charge 
qubit, dephasing processes are often neglected by biasing 
the qubit at the degeneracy point. However, as shown 
in Fig. [31 at the optimal bias point (sq s» 1 GHz) , the 
pure dephasing rate — 8 MHz is larger than the re- 
laxation rate = 4 MHz. Hence, it is important to 
include dephasing in the study of ground-state cooling. 
Ground-state cooling of a MR can only be realized if the 
qubit dephasing rate is sufficiently low, see Fig. [3] 

In conclusion, we have discussed the quantum theory 
of the cooling of a MR using the back-action of a super- 
conducting flux qubit. We have shown that the present 
noise level of the qubit allows ground-state cooling of 
the MR. The approach used in this paper can be ap- 
plied to other qubit-assisted back-action cooling schemes 
to estimate the cooling limit under the influence of pure 
dephasing. Our method can also be applied to different 
qubit noise spectra. 
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